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ABSTRACT 


A mathematical model was developed to include line 
effects in the submarine hydraulic system dynamic 
performance analysis. The project was undertaken in an 
effort to demonstrate the necessity of coupling the entire 
hydraulic power network for an accurate analysis of any of 
the subsystems rather than the current practice of treating 
a component loop as an isolated system. It was intended 
that the line model could be coupled with individual elements 
in modular computer programs to facilitate design and alter- 
ation decisions. 

The line model was developed and coupled with a submarine 
driving-system servo-loop in order to demonstrate the above 
premises. However, due to programming difficulties 
additional work is necessary to successfully complete the 


analysis. 
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I. INTRODUCTION 


A. GENERAL DISCUSSION 
The vast majority of hydraulic systems currently in use 
are of a dedicated type. That is, they are designed for a 
specific function with consistent and well defined constraints. 
In most cases they are further limited to a loop in which 
the pump supplies power to a single power element. Such 
systems can be delivered as compact units with a minimum of 
connecting piping. Some examples of dedicated systems in 
Naval use are: surface ship steering units, gun mount 
train and elevation units, and ammunition conveyor systems. 
Although submarines have employed dedicated systems in 
the past, all moderm U.S. submarines use constant pressure 
systems [Ref. 1, Tab. 561-1]. The constant pressure systems 
contrast with the dedicated type in that the supply and 
return headers extend the entire length of the ship and serve 
a multitude of power elements. The total piping length of 
these systems may exceed the length of the ship by a factor 


of four or more. 


B. THE NEED FOR A LINE MODEL 

For the packaged configuration of the dedicated system, 
with its relatively short piping lengths, the current practice 
of neglecting line dynamics and lumping the system volume 
has produced acceptable analytical results. However, even 
in systems with low line-volume to total-volume ratios, the 


computed values of the hydraulic natural frequency (wh) 


Ll: 








are frequently in excess of 40% higher than measured values 


[Ref. 2]. Since w, sets the upper limit on system response 
and is affected by the volume as well as line dynamics, a 
model which includes these factors would undoubtedly produce 
more accurate results. 

It was the intent of this project to demonstrate that 
the submarine hydraulic system, with its extremely long lines, 
requires that the line dynamics be included to achieve any 


degree of accuracy in analytical solutions. 


Cee POTENTIAL FOR SYSTEM MODEL 

Discussions between the author and Mare Island Naval 
Shipyard Hydraulic Design Division personnel (code 260) , 
as well as a Naval Sea Systems Command representative 
[Refs. 2 and 3], disclosed that there is currently no computer 
program to model the dynamic behavior of a submarine hydraulic 
system. Such a program would be useful in the design of new 
systems, both for sizing of lines and elements and for 
determining the response of components served. For existing 
systems, it would allow a determination of the effects of 
alterations prior to their implementation. 

Additionally, a computer program of a dynamic system 
model would be advantageous in analyzing vibrations trans- 
mitted through the system. The ''noise' generated by 
vibrations is a critical element of the "Submarine Sound 
Quieting Program,'' while the fatigue due to vibrations is 


important to the life of mechanical elements in the system, 


especially the seals [Refs. 5,6, and 7]. 





The development of a dynamic line model which could be 
programmed in a modular format and coupled with models for 
the individual elements would facilitate the creation of 


such an overall system model. 


i 





It. SYSTEM MODELED 


A. MAGNITUDE 

It was recognized early in the undertaking that it would 
not be possible to develop a model for an entire submarine 
hydraulic system in the time available. For this reason, 
the stern diving planes servo-loop was chosen for the 
primary development (Fig. 1). This sub-system was considered 
to be representative of one of the more complex loops in the 
overall system. ( a loop being defined as the fluid path 
through a control valve - to the power element - back 
through the control valve). Successful modeling of the 
stern planes servo-system would readily permit expansion 
of the system to include the fairwater planes and rudder 
loops in parallel operations. With these three sub-systems 
modeled, the remainder of the hydraulic system could be 


relatively easily developed. 


bee REALISM 

An operational ship's hydraulic system was deemed 
essential as the basis for the model, as it would allow the 
computed parameters to be compared with actual performance 
data. The hydraulic plant of USS PERMIT (SSN 594) was 
chosen as the system to be modeled. It was considered 
representative of a modern submarine hydraulic system and 
the technical information was readily available at Mare 


Island Naval Shipyard. 


14 











WALSAS SANVTId NYALS 


GVaHSSOUO 


9 a = 





f JN ITANOD 
<— X oe 
SUNVId NWALS 
TINH 
jTuaSsaMd 


A’IlddNsS 






SANT TAO 
yy WVd 





TUTTI TY 


om 


NaN Lad 


—OLINVUGAH 





AOVE -Cadda 





T ean3sty 


NOILVIS 
"LOULNOO 





f 
wey \ H 
YN 


15 





ieee 1 DUO UL D 


Bee LYPE 

The hydraulic oil used in all U.S. submarine internal 
high pressure (3000 psi) systems is Navy symbol 2190TEP. 
This is a petroleum based fluid conforming to military 
specification MIL-L-17331 [Ref. 8]. It exhibits 
approximately the same characteristics as an SAE-30-Detergent 


motor oil. 


foe PROPERTIES 
1. Density 
The density (9) of the oil was obtained from Ref. 9 
in terms of specific gravity (SG). The values supplied 
were 
0.896 at 74°F. 
0.0901 at 60°F. 


SG 


SG 


The density of water at 60 degrees Fahrenheit is 
pea /5 x 107? (slug/in?) and is assumed to be the reference 
for specific gravity unless otherwise specified [Ref. 10]. 
Specific gravity is defined as the ration of the density of 
the oil at the given temperature to that of water at 60 degrees 


Fahrenheit. 


SG = p/p (1) 
0° "Wed 


Solving for the density of oil at temperature (T) 


16 





yields 


S422 i To? slug/in> at 60°F 
5 


p 


oo 8.3754 x10 slug/in> at 74°F 


2. Cubical Expansion Coefficient 

The cubical expansion coefficient (a) is the fractional 
change in volume of the fluid due to a change in temperature and 
is used in the determination of density, viscosity (u,v), 
and bulk modulus (8) values. 

Density values arrived at in the previous section 
were used to obtain a value for the cubical expansion coefficient. 

Using the form given in Ref. ll, 

e(P,T) = o [1 + (P-P,)/8 - a(T-T,)] (2) 
where the subscript o) denotes reference values of one 
atmosphere of pressure and 60°F temperature. 

Since the values for the density of the oil were 


obtained at one atmosphere 


P = ie 
equation (2) reduces to 
peewee ite a (T - T)] (2a) 


Inserting appropriate values into (2a) and solving for a 


yields 


ome 


a = 3.9638 x 10° OF 


me Viscosity 


a. Kinematic Viscosity 
Reference 12 specifies that 2190TEP oil shall 
have a kinematic viscosity (v) of 190SSU (Saybolt Seconds 


Universal) at a temperature of 130°F. 


ky 





Using the standard conversion [Ref. 14} for 
kinematic viscosity, where the time (t) is greater than 
100 seconds (SSU) 

OP 2Z20REe= 4857 t) (3) 


Vv 


yields 5 
41.09 cs at 130°F 


) 
v = 92.08 cs at 100°F 
Converting from centistokes (cs) to English units was 
accomplished with the conversion factor (Ref. 15) of 1.549 *107° 
(in«/sec)/cs. Thus the values obtained for kinematic viscosity 
were 
v = 0.06367 in*/sec at 130°F 


v = 9.1427 anaes at 199°F 


b. Absolute Viscosity 
Knowing values of density and kinematic viscosity 
permits the determination of absolute viscosity (y). 
Hw = pv (4) 
Using the form presented LTS emer. 


Meee - A (Tt - 7.) (5) 


O 
at one atmosphere of pressure, where A is the temperature 
coeffcient of absolute viscosity with units of (OF ty. 
Combining equations (4) and (5) and then 
solving for lambda gives 
ee amo py) Gar) (6) 
Arbitrarily setting 
T, = 100°F and T = 130°F 


in conjunction with values of density and kinematic 


viscosity at these temperatures yields 


18 





T, = 100 T = 130°F 
p, = 8.2886 x10? sp = 8.1884 x10 
vo = 1.427 x107+ seeeG eo Gna 0as 


Substituting these values into equation(6) produced a 
value of 

\ = 2.731 x107* OF t 

Reference 17 presents the pressure conversion 
for absolute viscosity as 

logjo(u/ul) = cP (7) 
However, as the other functions are related by the difference 
in pressure divided by the bulk modulus (AP/8), it would 
seem reasonable to assume that c in equation (7) has a form 
of (c,/8). This would then yield an equation for absolute 
viscosity at constant temperature of 

H=ule -(c, 4P/8) (8) 
Lacking any other information, Cc, was set equal to one, 


and a final form for absolute viscosity was formulated as 


u = ue((P/8) - A(T-T,)) (9) 
where 

Bas Poe. 

P = reference pressure (o psig) 

P = system pressure (psig) 

8 = system effective bulk modulus (psi) 

T = system temperature (°F) 


T.= reference temperature (°F) 


ee 





4. Bulk Modulus 
Although the literature on hydraulic systems indicated 
that the bulk modulus (8) would be supplied with the oil, 
visual inspection of shipping containers and bills of 
lading proved futile. Likewise, the military specifications 
(MIL-L-17331) failed to produce any worthwhile information 
in this regard. A nominal value of bulk modulus for the 


oil was therefore assumed to be 220,000 psi as recommended 


3 


mere r. 18. 


C. LIMITATIONS 

In a further effort to maintain the scope of the project 
within reasonable constraints, it was decided to treat the 
temperature as constant. Inferences drawn from the Naval 
Ships." Technical Manuals previously cited and the author's 
own experience, indicated that a temperature of 130°F was 
a reasonable choice. With this value of temperature, the 
equations for the other properties were simplified. 

From equation (2) 

9 = 8.188 xl0~> + 8.4221 x107>(P/8) (10) 
and equation (9) yields 


Me 513 x10 ° 5 (P/8) . (11) 


Zo 





IV. » MODEL DEVELOPMENT 


A. SYSTEM OPERATION 

Submarine steering and diving control surfaces (rudder 
and diving planes) are positioned through mechanical linkages 
by hydraulic rams. Each surface has normal as well as 
alternate sources of hydraulic and electrical power to ensure 
reliability. For the purpose of model development, only the 
normal mode of the stern planes was considered. 

The normal mode of operation for the SSN 594 class 
of submarine, employs a position control electro-hydraulic 
servo-loop (Fig. 1). In this type of system (Fig.2), control 
station operator manually positions a stick to achieve a 
desired angle of inclination of the planes. This produces 
an electrical signal which is summed with the feedback 
Signal from the position output to create an error signal. 
The error signal is amplified and transmitted to the torque 
of the servovalve. The torque motor displaces the wand spring 
of the servovalve-pilot, which causes an imbalance in pilot 
pressure, which in turn shifts the main spool allowing the 
servovalve to pass oil to position the ram. As the desired 
position is attained, the error is reduced and the pilot 
returns the spool to a neutral positions, blocking the flow 


Gn O1l. 
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B. SERVO-SYSTEM 
The area of interest for the model development was the 
hydraulic fluid flow in the servo-loop as depicted in Figure 3. 
By inputting the valve spool motion and ram load force as 
explicitly defined functions, the influences of the external 
factors involved in the feedback loops (the servovalve torque- 
motor and pilot valve, and the ram loading) were minimized. 
This permitted the generation of flow and pressure equations 
around the loop with a minimum of complications due to the 
ancillary characteristics. 
1. Servovalve 

The two servovalves commonly in use in the Submarine 
Force are both manufactured by the Sargent Corporation and 
are designated Sanders Model SV 438-10P and SV 438-15P 
[Ref. 19]. Design and performance of both valves are 
Similar. The SV 438-1LOP which is installed in the (SSN-594), 
was used for this simulation. Figure 4 is a representative 
drawing of the actual valve employed - a four-land four-port 
element. However, comparison with the valve depicted in 
Figure 3 - a three-land four-port model - will show that 
the flow is essentially the same in both valves. It was 
considered easier to follow the flow paths of the valve in 
Figure 3. 

All technical data for the servovalve was obtained 
Peom Ref. 20. 

a. Flow Paths 

The pilot valve flow was disregarded,as it has 


its own supply and return lines. 


a 





RETURN 


SERVOVALVE 


5 iN 74 aa 
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Designating the supply flow as Qin? the return 
flow as Q-» and the other flows by their respective nodal 
points (Fig.3), the flow equations were derived. Using 


the form given in Ref. 21, 


BO eee = dv,/dt sts (V,/8) dP/dt (GiZ)) 
the results were: 

ry ea gs ) 

Oe eo tt Oy, (14) 

ee armas ae, / 6, + V, (15) 

QQ = & + % + VoPo/By + Vo. (16) 


Here the volume across the lands of the valve was approximated 
as zero. The superscript (°) represents differentiation 
with respect to time. 
Since the valve is manufactured from high 
strength alloys, 
i oo 


V, = Vp, 


iE 
as shown in Appendix A. 

In addition to summing the flows through the 
valve, the flow rate was also defined through the lands for 
three different valve positions. 

(1) Null Deeimions. If the centered (null) 


position, a critical-center servovalve will have leakage 


past the lands.’ This flow is described in terms of null 


‘For a comprehensive coverage of all terms used in this 
paper, the reader is referred to Ref. 2. 
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flow pressure coefficient an as 


Q = 2K CAP). (17) 
The four node flows for Xx, = 0 were 

Q) = Ko (Pi y-P) (18) 

Qo = BK 9 (Pin Po) (19) 

Q3 = ¥K,,(P-P_) (20) 

Q, = BK Og (P27 PL) C21.) 
with 

P., = Py = Pe 


See appendix B for values ie used for the model. 

(2) Valve Spool Off-Centered Left. The valve 
travel (X) was arbitrarily defined as positive to the left 
(Fig.3). Employing the orifice flow equation [Ref. 22 } for 
the fluid passage through nodes one and four, and an empirical 
relationship for nodes two and three, the flow equations were 


determined to be 


Q) = CgAY 20Pi "Py d/o e-) 
Q = BK e-A(K,/X 0)” (Ppa Po) (23) 
Q, = RK e-A(X,/K)° (P,-P,) (24) 
Q, = CgAY2(P9-P,)/P, (25) 


where om = 0.61 for full periphery rectangular ports and 

A = md, XX, (26) 
with dy defined as the land diameter of the spool valve. 

Since a value of leakage flow was not 
determinable from Ref. 20, the empirical relationships of 


Appendices C and D were developed using Ref. 23 as a guide. 


‘sf 





(3) Valve Spool Off-Centered Right. The 


equations for the "negative" flow direction were similar to 
those developed in the previous section, with the obvious 
exception that Qo and Q3 became the principal flows. In 
order to maintain their flows as positive values, Xx. was 
squared and included in the radical. 


The form for Xx, less than zero was 


Qy = Bape M(Ky/ Xm) (Pin -Py) (27) 


Q = Cg A 2K (Ps a-Pp)/0y_ — 

Qg = Cguv! 2K (P]-P,)/09 (29) 

Q, = #K,(e-A(K,/X,.) (Py-P_) (30) 
with 

Wo = nd, (31) 


and all other variables as defined in previous sections. 
2. Lines 

The line equations were generated with the intention 
of accounting for all factors influencing the pressure and 
flow in the system. The schematic diagram of the actual 
system employed in this development is shown in Figure 5. 
For development purposes, only the odd numbered nodes are 
shown here as the even numbered side equations were identical. 

a. Flow And Pressure 

Disregarding all but the main line piping, the 

flow path from the servovalve to the ram was as depicted in 
Figures 3 and 5 were identified by the same numbers. The 
dynamics of these segments were then defined at each point 
when viewed as in Fig. 6. 
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FIGURE 6 LINE NODE IDENTIFICATION 


Point A corresponds to nodal point (5) of Figure 3 
and is considered to be located at the entrance of the pipe 
at its junction with the servovalve. Similarly, point B is 
located at the connection of the two different diameter pipes, 
in the entrance of the downstream pipe. It is equivalent 
to node (7) of Figure 3. Location C is just beyond the end 
of the line segment at the entrance to the cylinder, 
identified as node (9) in Figure 3. 
Utilizing these definitions for the nodal points 
and viewing the direction of flow as always being from A to 
C, the dynamic equations for each element were developed. 
(1) Node Five. The pressure and flow equations 
for node five were derived from Ref. 24. 
Qs = O7 + V57P/Bs7 + V5, rs? 
Pe = P, - pcKaQs/2Az (33) 
The double-subscripted variables were employed 
to denote average values through the line segment, while the 
Single-subscripted variables represent values at the particular 


node. Thus 
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V = AL-L 


SVE Geb 57 (34) 
with Lo5 being the actual length of pipe from node (5) 
to node (7) and 


hE mDS/4. (35) 


5 
As can be seen in Figure (6), the diameter of segment A-B 
(5-7) is constant and equal to the diameter at (5). The 
total volume and actual length will be developed further 
in later sections of this report. 

The average pressure in the interval was 

esp ea Po)/2 (36) 


with , 
Poo = JP.7dt . (7) 


The change in volume with respect to time (V) 
was set equal to zero for all line segments, because the 
lines are made from high strength materials which do not 
permit any significant expansion or contraction. The 
exception to this would be in cases where there are flexible 
sections installed ( Aeroquip -type hoses). There are 
numerous locations in the system where these are installed, 
but not in the servo-loop modeled. In the instances where 
flexible sections are incorporated in the system, an expression 
for the change in volume would have to be developed based on 
the properties of the material used. 

The value of Ke was based on an abrupt change 
im the cross-sectional area of a pipe,and the development 
is shown in Appendix E. As used in the above equations, the 
Sign of Ke changes with the direction of flow such that the 


second term of equation (33) will always reflect a pressure 
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loss for the referenced point. Thus, when flow is positive, 
as shown in Figure (6), Py will be greater than Pe, but when 
the flow is reversed (negative in Figure (6)) Pe will be 
greater than P,. 

Qz is the flow at node (7) and will be defined 
in a subsequent section. Pi 1s the pressure in the main 
servovalve cavity and is 


P, = {P dt (38) 


where Py was introduced in section IV.B.l.a. 

(2) Node Seven. The flow through node seven is 
exactly the same as that through node five, with subscripts 
changed to reflect the appropriate line segment. 

7 = Q% + VigP7q/Brq + V7q 7 

Pog = JProdt. 

The pressure at point (7), on the other hand, 
must reflect the losses due to flow in line segment (57). 
This was accomplished with a combination of the equations 
in Ref. 24 with those in Appendix E: (40) 
Py = Ps - 32ug7bes7Q57/AsD5 - 1.14957057/A5 -07K7Q7/2A5 

The second term on the right side of equation 
(40) is the Hagen-Poiseuille term for fully developed laminar 
flow in pipes. The third term accounts for developing laminar 
flow. The final term accounts for the entry/exit loss at 
cross-sectional area changes, as shown in Appendix E. 

All three terms must reflect a loss of 
pressure for flow in either direction. Since oer assumes 


a negative value for flow reversal, it presents no problem. 
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However, ae, will always be a positive quantity unless it 
is defined as 

Qe = Geo [Qcal, (41) 
in which case the sign of One will reflect the direction 
of flow. As noted previously, Ky determines the sign of 
the final term. 

Az is used in the last term of equation 
(40), becaue the area employed is always that of the smaller 
of the two at a junction. 

Previously undefined variables are 

Qs7 = (Qs + Q,)/2 (42) 

Les, = effective length from (5) to (7). 
Effective line length will be discussed in a later section. 

(3) Node Nine. The pressure at point nine has 

the same form as that at seven 
Pg = Py-32uzglezgQ7q/AzD7 - 1-14079Q59/A7-PgK9Q5/2A7 oar 
with 

Q59 = Arq 1059) (44) 
for the reasons discussed in the previous section. 

Flow through this junction is 

89 = [11 F er, ey 
Q,, will be defined with respect to ram flow in a following 
segment of the paper. 

Soeme oto) (46) 
The coefficient of external leakage (C,1) is the allowable 
leak rate past the cylinder tail rods to the atmosphere. 
Appendix F derives the values used, based on information 
contained in Ref. 25. 
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There were no volume terms included in 
equation (45) because the transition from the end of the 
line to the cylinder entrance was considered an infinitesimal 
distance. 
b. Secondary Component Effects 

It is readily apparent from a comparison of 
Figure (5) with Figure (3) that the simplified system omits 
numerous valves, fittings, and other miscellaneous elements 
which are installed between the servovalve and ram. Although 
these "secondary ‘components have a minor influence individually, 
when combined they contribute significantly to both static 
and dynamic characteristics of the system. 

The approach taken in this development was to 
consider them in three categories: branch lines, minor 
components, and major elements. 

(1) Branch Lines. Branch lines are the segments 
of pipe connected to the main line. They are subjected to 
system pressure but do not ordinarily allow fluid to flow. 
Gage lines, hand positioning branches, and capped sections 
were included in Figure (5). 

These lines contribute only to the total 
system volume. 

(2) Minor Components. The minor components 
were considered to be any elements installed in the system 
which contributed to the effective length and volume of 
the system but were not specifically listed in the list of 


major components [Ref. 26]. 
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The elements in this group were ball valves, 
tees, and elbows. Reference 27 was used as the source for 
information on them and Apvpendix G explains their development 
for use in the overall analysis of the system. 

(3) Major Components. The deceleration valves 
(decel) and the blocking valve were included as major components. 
These were listed in Ref. 26 with specific pressure and flow 
characteristics. The volumes employed were a "best estimate” 
based on diameter and the author's personal experience. The 
effective lengths were as determined from the procedure in 
Appendix H. 

The function of the blocking valve is to 
transfer to an alternate source of hydraulic supply (the 
Vital System), in the event of a loss of main header pressure. 
As pressures in the two systems are ordinarily the same 
(within 1lOOpsi. of each other) and the valve has internal 
seals, cross-port leakage is negligible. 

The deceleration valve is a cam operated 
full-flow valve which is installed to prevent over-travel 
of the ram. As the ram approaches its maximum limit of 
travel, a collar on the tail rod mechanically trips the 
decel can, shutting the deceleration valve. This permits a 
fluid cushion to arrest the ram travel rather than allowing 
it to a "hard-stop" (metal - to - metal) at maximum velocity. 
Motion in the reverse direction is permitted by flow through 
the internal check valve until the linkage allows the main 


deceleration valve to reopen. 
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(4) Combined Effects. The various elements 
were combined in tabular form for inclusion in the system 
analysis (Tables 1 and 2). These were input directly into 
the computer program to solve for the various parameters 
as discussed above. 

3. Cylinder and Piston 
The final element in the system is the cylinder and 
piston combination. Here the power is transmitted from the 
fluid to the stern planes by the piston and its associated 
linkages. 
a. Fluid Dynamics 

The pressure and flow equations assumed a form 

Similar to those of the lines. 


(1) Flow. The flow in the cylinder was described 


by 

Ce ee a aa! 87 (47) 
where 

Qag = Sp e(Pyy7P 12) (48) 
with . 

age JP, dt (49) 


Q; 2 is the internal leakage from the high 
pressure side to the low pressure side of the piston. The 
determination of Ci p is shown in Appendix F. Flow was 
defined as positive from node (11) to node (12); thus,when 
Pio is greater than Paw Qs will assume a negative value 


as required by the line equations. 
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Aqy is the effective piston area 
Z 2 ; 
where Di4 is the land diameter of the cylinder and D.. is 


the diameter of the piston tail-rod. 


r is the velocity of the ram. 


Yy1 = Arita (52) 
where 

Ly, = 4, +X (93) 
Pad Ly = (Ly -tp) (54) 


Ley is the length of the cylinder Cy 
is the thickness of the piston. Ram position was defined 
meepositive to the left in Figure (5), so that motion to the 
left increases Vii while movement of the piston to the right 
decreases Vi1° For this reason, it was necessary to define 
Vio = L 2 (55) 
and ultimately 


Qi2 = -Q9ie - VYi2 - Vu2Pi2/8i2- (36) 
Thus, the various flows in the system would automatically 
assume the proper sign to reflect the direction of flow. 

(2) Pressure. The pressure equation for the 
cylinder and piston arrangement became 
- ji Z Z 

Par ~ Po ~ 32¥yyby 3% /4ppPyy, ~ 2-14 9471/47) - (97) 

There was no entry loss term because it 
was taken into account with Pg- The other difference, with 


regard to the previous line loss equations, is that Dz4is the 
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@eanad diameter of the cylinder while A, , is the effective 
area of the piston. The effective area was used because 
Q/A is a velocity term employed in defining viscous loss. 
On the other hand, Di was used to account for the surface 
area upon which the fluid is "dragging"; i.e., the walls of 
the cylinder. Since the ali rod is traveling at essentially 
the same velocity as the fluid, it contributes no losses 
due to viscosity. 

As in previous sections, it was necessary 


to define 


OF) = 7°11! (58) 
to account for flow direction. 
b. Piston Dynamics 
Ram motion was defined in terms of position (XL), 


velocity Cor and acceleration (X_). 


ee part (59) 
k = fide (60) 
X= (Az Py 1-Ap Py 2-5 phe Fy) /Mp (61) 


The load force (FL), as described at the outset, 
was input as a ramp function dependent on ram position and 
velocity. 

The eller of oil (M34) term was included to 
account for the force required to maintain the momentum of 
the oil. It was determined to be 

141° 012%19 SOO (62) 


and . 
for x <0. ; (63) 


Mo¢1 7911"11 
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The total mass of the system (Mp) was defined 
as the mass of the ram plus the mass of the oil from the 


low pressure side of the piston to the return tank, such that 


tank 
M,. =m + y Cae (64) 


ve ey 
ram 
where the subscript (;) indicates each main line segment 


(flow path) from the ram face to the return tank. 


C. PROGRAM 
The foregoing equations were re-arranged to yield the 


form shown in the attached computer program. They were 


modeled in a Continuous System Modeling Program (CSMP IIT) 


in the belief that it was the best modeling method available. 
However, numerous attempts failed to accomplish the desired 


mesult. 


The algebraic equations form a loop for which no solution 


was found. 
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V. SUMMARY 


A. CONCLUSIONS 

It is felt that the development of the line equations 
demonstrates the necessity of including these factors in 
any in-depth analysis of the submarine hydraulic system. 
It is also the author's belief that the approach taken in this 
presentation will yield worthwhile results for any system 
employing extensive quantities of piping. 

As stated in the introduction, a generalized computer 
program of this type would have many benefical aspects for 


both the submarine designer as well as the operator. 


B. RECOMMENDATIONS 

The first recommendation is a possible means of clearing 
the algebraic loop problem currently preventing the program 
from executing. This would require clearing all of the 
pressure and flow equations of reference to any variables 
except for themselves or ones which directly invoke an 
integral. That 1s 

Po = £(P59,P55,P,,P5,etc.) . 
However, it could not take the form where Po was set equal 
to Pg or any of the flow terms (Q; ), as they do not invoke 
an integral function. 

The other possibility would be to locate another program 


which could efficiently handle both integral as well as 


Simultaneous algebraic equations. 
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APPENDIX A 


Determination of Servovalve Volume 


No values were obtainable for the volume of the 
servovalve. Therefore, the following method was employed 
to obtain a realistic approximation. 

Using nominal dimensions given in the technical manual 
(Ref. 20 Jand the knowledge that the valve is made of 
aluminum alloy, a volume to weight ratio was determined. 

Dimensions ion) x G40025ex99.5 in. 
Specific Gravity of alloy approximately 2.8 
Given weight 69 lb. 

Density of water at 60°F. 62.41b/£t> 


_ cee 
Voody = 1417 in 
Yalloy = 0.101 1b/in” 
W = Vy = 143.3 lb for a solid block 
then 
Weiven! = 0.48 . 


Thus the valve body weighs approximately half of what a 
solid block would weigh. It was therefore assumed that 
the volume of metal in the valve was 708in?. Further, 
estimating that the passages for the pilot system occupy 
half of the void area yielded a main valve volume of about 
350 in>, and since this volume is that occupied by both the 


supply and the return segments of the valve, half was assigned 


to each. 





APPENDIX B 


Null Flow Pressure Coefficient Determination 
The allowable null flow leakage was given as 600 cc/min. 
at 1200 psi. [Ref. 20]. This was converted to standard form 


of (in?/sec)/psi. as follows 


K.. = 600(cc.min)/1200 psi x 1 in?/16.39cc x 1 min/60 sec. 


This yielded 
Ko = 5-0844 xl0* (in®/sec) /psi. 
and 


aK. = 2-5422 xl07* (in?/sec)/psi . 


APPENDIX C 


Empirical Valve Leakage Coefficient 


The leakage flow across a valve port decreases as 
the valve moves to block the port. Since a plot of leakage 
flow (Q,) versus valve position (X.) is symmetrical, with an 
exponential decay as shown in Fig. C-l, and exact data was 
not given for leakage other than the null position an 
empirical expression was developed. 


Qe 
Q, MAX 


Figure C-l Leakage Flow vs Valve Position 
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The most appropriate expression to fit this curve took 


the form of 


_ Z 
oe e7 + (Xy/ Ka) 


where 


Q, = leakage flow at any valve position 


— = L 
COA maximum leak rate aK Ao 


X maximum valve travel from null. 


Vmax 
Representative values were then chosen such that 
pee eat XO = 425 
That is, at 25% of maximum valve, travel leakage through 
the "blocked" port of the valve decreases to 2% of the null 
position leakage. 
Solving for A 
h = (1N)Qp ray! Qp))/ (Ky/ Xyrnaase) 


A = 62.6, 


2 


In reality A could vary from near zero to about 100. 
Realistically it would have a value between 50 and 100. 


4 = 65 was chosen as the value to employ in this 


presentation. 
The form used was 
IE 
Qp = aK ey! X ymax? 


with the determination of maximum valve travel (XD shown 


in Appendix D. 
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APPENDIX D 


Determination of Maximum Valve Travel 


Reference shows that 


Qrae = 80 GPM = 80 x 231/60 in?/sec 


at 
T = 115°F 


AP = 1200 psi 


total 


A = 600 psi 


Poort 
as depicted in Figure D-1. 


0 i PeQum rol 


E 


a 


po” Ss 
600 PSI 


BuetieewU=) Valve Tes, Flow Path 
Using the orifice equation 
3 
aS C A(24P/p)* 


with 


OQ 
ll 


amen 
8.442 x107? 


A =nd,X,, where dy =—s 1 in. 
X= X 


V vm 


Q = 
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such that 
Koen = (Qq/ mC gdy (24P,,/0)%) 
resulting in Kal 0.0426 in. 
This value was used in spite of the fact that it seemed 


large for a “standard" pilot controlled servovalve. 


—————— 





APPENDIX E 


Loss Coefficient Due To Abrupt Area Changes 


The standard form of the pressure loss due to abrupt 
piping diameter change is given as 


H, = KV"/29 


il 
where 
Hy is head loss, expressed in feet of fluid 
Weis Prutd velocrty 
g is acceleration due to gravity - 
Head loss can also be expressed as 
eS AP/a 
Y = 0g 
and 
y 1s weight density 
while velocity is 
V= Q/A = 4Q/nD° 
combining the above equations 
AP = SeKOovn“D" . 
K is defined from Figure D-1l 
*(1-D,/Do) for +Q 
-(1-(D,/D,)*)* for -Q. 


K 


K 


NOTE D, is always the smaller of the two diameters being 


considered. ee Ct 
De Dy 


Figure E-1 Flow at an Abrupt Area Change 
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APPENDIX F 


Seal Leakage Coefficients 


Seal leakage was based on permissable values listed in 
Naval Ship's Technical Manual,[Ref. 25.] The values used in 
this presentation were for "used" seals with maximum allowable 
leakage before replacement is required. This was felt to 
be reasonable since the operating temperature of the oil 
is normally greater than the test temperature, which will 
produce slightly higher lead rates. 

The leak rates were given as one milliliter per ten 
cycles per inch of rod diameter for external leakage, and 
five milliliter per five minutes per inch of rod diameter 
for internal leakage. These were based on normal operating 
pressures and a fluid temperature of LOO°F. 

Operating pressures were assumed to be 


P 


a 2000 psi 


P 1000 psi 


ae! 
where Pal 1s the pressure for external leakage and P., is 
the pressure for internal leakage. A cycle was defined 

as the time to go from neutral (zero plane angle) to full 
rise to full dive and back to neutral position. For the 
system modeled, stern planes travel is plus or minus Den 
Thus,a full cycle would be four times this value, or 108°. 
The normal rate of travel is about five degrees per second. 


Thus, 1108/5 = 21.6 sec/cycle or approximately 


3 cycles/min. 


oe 





Since the tail rod used, had a diameter of four inches, 


the allowable extermal leakage was four milliliters per ten 


cycles. 
Then 
aml , lec y SEY = 1, 2cc/min 
lOcy Iml 1 min 
and 
Cop = 1.2 cc yl min, l in? se i 
min ~ 60 sec .39cc 2000 psi 
aoe x107/ (in?/sec)/psi 


Similarly the internal leakage for a land diameter of 
8.25 inches, rounded to nine inches, was nine cubic 
centimeters per minute. 


And 
» 9 
Ci, * 60 x16.39 x 1000 


es oD sao (in?/sec)/psi 
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APPENDIX G 


Minor Component Analysis 


The values used in the derivation of the minor component 
effects, considered herein, were obtained from Crane Technical 
Paper No. 410, [ Ref. 27]. 

The volume was determined to be 

V = AL 
where A was the cross-sectional area of the pipe entering 
or leaving the element and L was a value determined for 
each type of element. The values for L were chosen as 

Ball valves te= 3 x D 

Tees L=2xD 


Elbows L=1.6 xD 
In cases where the element has differing entry and exit 
diameters, D is the larger of the two. 


The values for the effective lengths were 


Ball valves ee 3 x D 
Flow through tees L, = ZO x D 
Branch flow tees Ee = 60 x D 
45° elbows L, = Oe ) 
90° elbows L, = 30 x D 


where the diameters were as above. An additional factor of 
K (appendix E) would have to be included for diameter 
differences in the main flow path. 

The value of effective length cited above, were determined 


for Reynolds numbers (R,) greater than 1000. Therefore, 
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it was necessary to modify them 


Since 
and 


with 


then 


Et 


So that for 


Or 


and for 


R, = 4Q/ (Dv) 
Log L,.R,/1000 Lor, R,<1000 
Log the desired effective length 
Luke the calculated effective length from above 
= -3 
ie = be (4Q/(1rDv) ) x10 
R.>1000 then L =L 
e- é ec 
4Q/mDv>1000 
Qo/Du> 785.4 t= Joe 
Qo/Du<785.4 L, = 1.273 x107°L, .Q0/Du 
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APPENDIX H 


Effective Length of Major Valves 


The standard head loss equation 


oo = Ko /2g was rearranged in the form 


L 
AP/pg = KQ*/2gA° 
This yielded 
K = NDE 
since 
A =nD*/4 
K = (APD? / 800° 
and since 
a L/D 
L. =1-APD?/8 90° 
To allow Q to be entered in the form presented (GPM), 
a conversion was included. 
231 in>/gal x 1 min/ 60 sec 
so that 


Q (GPM) x 231/60 = Q (in-/sec) 


and L, = .0832 PD? /Q° 

where 
D = diameter of valve (in) 
P = given pressure loss (psi) 
Q = given flow (GPM) 


density of fluid at given temperature (slug/in?) 


ie: 
Il 





APPENDIX I 


TO EMPLOY: 
SUPPLY LINES; Dl IS ALWAYS THE POINT CLOSEST TO THE 


PUMP - THE UPSTREAM POINT 


RETURN LINES; Dl IS ALWAYS THE POINT FARTHEST FROM 


THE TANK - THE UPSTREAM POINT 


DUAL DIRECTION (REVERSING) LINES; Dl IS THE POINT 


ENDMAC 
INITIAL 


CONSTANT 


CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 


CONSTANT 
CONSTANT 


CONSTANT 


CONSTANT 
CONSTANT 


CONSTANT 


CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 


NEAREST THE CONTROL ELEMENT (VALVE), VIEWED AS IF 
THE LINE ALWAYS SUPPLIES PRESSURE TO THE POWER 
ELEMENT 


KcCccc=600.,FT=1200. , LAMB=65. 


VALVE LEAK RATE OF MAIN SPOOL 
NOTE: LAMB IS STRICTLY EMPIRICAL AND BASED ON AN 
AVERAGE CURVE OF LEAK-RATE VS VALVE DISPLACEMENT. 


XVM=0 .04.25, DLV=1.0, CD=0.61, V1=175. , V2=175. 


ODS=1.90, THS5=0.25, L5=l.8.0, NTB5=0. ,NTT5S=0.;5... 
NLU5=0.,NL95=1. ,NV5=0. 

po eos to —0,25,L0=10.0, NTBO=-0. ,NTTO=0.,... 
NLY6=0.,NL96=1. ,NV6=0. 

Siya2 0, by -O2e7, L(=cee0.0,NTB7=c. ,NTT7=1. , NLU7=0.,... 
MEov—cee NY 7=1. 

op8=2.0, TH8=0.27, L8=216.0, NTB8=3. , NTT8=0. , NLIB=0.,... 
NL98=5.,NV8=1. 


OD71=0.75,TH71=0.109,L71=8.0, NTB71=1. , NTT71=C.,... 
ime l=0.,NL971=0.,,NV7L=0. 

Gome-O0 57 5,647 2-0-0050, e7e=le2O-0,NTH72=0. ,NTT72=0.,... 
NLL.72=0.,NL972=1. ,NV72=1. 

OD7 3=1.5, TH7 3=0.203,L7 3=12. ONTB7 3=0. , NTT73=0.,... 
NLL? 3=0. , NL97 3=0. , NV7 3=0. 


O@at—O. 7 Semued—0-1097081=6.0,NTBSi=1.,NTT81=0.,... 
Migtow—O se NELGOl=O.,NVGl=0. 

OD82=0. 3750, TH82=0.058, L82=120.0, NTB82=0. , NTT82=0.,... 
NLU82. , NL982= 2, ,Nve2=1. 

ObC3=1. 5, THS 3=0. 203, LS 3=12.0,NTBS3=0. , NTTS3=0.,... 
Nimiee—On0. NLC 3=0.,NV63=0. 


PDR7=60 .0,QDR7=50.0, RDR7=8. 68E-05, VDV7=60. , DDV7=1.09 
PBR7=50. ,QBR7=50. , RBR7=8. 68E-05, VBV7=20. , DBV7=1.6 
E=3.0E07, POIS=0. 3 

BO=200000. 

PR=50., PIN=3000., FA=0.1,O0DRE=2.0, THRE=0.18,VR-2750. 
2: O, CELCC=1.2, PCEL=2000. , PCIL=1000. 
DLR=8.255,DRR=.0, ODR=11. Ps LCY=39.875,TR=5.125,MR=12.1 
PI=3. 141593 
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KCc=0.5#KCOCC/ (PT#983.h) 
THE VALUE OF (0.5) IS USED BECAUSE THE TECH. MANUAL 
LISTS A LEAK-RATE FOR THE WHOLE VALVE; THE 
VALUE SOUGHT IS THAT FOR HALF THE VALVE (ONE SET 
OF PORTS). 


KCOCC IS TECH. MANUAL LEAK-RATE IN (CC/MIN). 
ay IS TEST PRESSURE USED TO ESTABLISH THE LEAK-RATE (PSI). 


KCO IS LEAK-RATE IN ( (IN ¢3/SEC)/PSI). 


AG=CD:: PI*DLV*SQRT (2. ) 


AG IS THE PRODUCT OF THE CONSTANT COEFFICIENTS USED IN 
IN THE SERVO-VALVE EQUATION. 

CFA=FA/1.) 

FA IS THE FRACTION OF AIR ENTRAPPED IN THE. SYSTEM, 


EXPRESSED AS A PERCENTAGE OF THE TOTAL SYSTEM VOLUME. 
CFA IS THE COEFFICIENT OF THE AIR TERM IN BULK MODULUS 
CALCULATIONS. 


K5P,K5M=ENTCO(D1,D5) 

K6P KOM=ENTCO (D2, D6) 
K7P, -K7M=ENTCO(D5, Dy) 
K8P ,K8M=ENTCO(D6, D8) 
KOP, KOM=ENTCO(D7,D11) 
K1OP,K1OM=ENT CO(D8,D12) 


THESE COEFFICIENTS ACCOUNT FOR EXIT LOSSES ONLY: 
ENTRY LOSSES ARE CONSIDERED AS EXIT LOSSES FOR 
THE NEXT ELEMENT Wee 


VP5,A5,LEP5, DSLAS, BCS5= VALE(PI, E,POIS,0D5,THSL5, NTTS,NTBS,... 
NLLS, NL95, NVS) 

VP6,A6,LEP6, D6LAG, BCS6=VALE(PI, E,POIS, OD, THO. Lo, NTL. NTBO,... 
NLL6, NL96, NV6) 

eA? ,»LEP,,D/,LA7 »,BCS/=VALE(P1,E,POIS,0D7,TH7 ,L7, NTT7, NTB7, o«- 
NLL » NL97, NV7) 

VP8,A8,LEP8,D8,LA8,B(S8=VALE(P1,E,POIS,0D8,TH8,L8,NTT8,NTBS,... 
NLUS ,NL9S, NV81) 


fA; BL, Le /BL,D/7B1,LA7B1,BCS7L=VALE(PI,E,POIS,OD71,TH71,... 
Pye Nite. NPS? L NLU? L,NLO7TL,NV7L) 
V7B2,A7B2,LE7B2 ,D7B2,LA7B2,BCS72=VALE( PI ,E,POIS,OD72,TH72, -.. 
L72,NTT72, NTB72, NLU72,NL972, NV72) 
V7B3,A7B3, LE7B3,D7B3, iA beyeecme—VAlbeorl, §,PO1S,0D/3,THi{3,-<. 
L7 3, NTT7 3, NTB7 3, NLU7 3, NL97 3, NV7 3) 
v8B1,A8B1, LESB1, D8B1,LA8Bl, BCS81=VALE(PI, E,POIS, Op81 PHO ligne lee 
L81, NP T81,NTBO1,NLUE1, NL961, NV821) 
V8B2,A8B2, LESB2, DBB2,LA8B2, BCS62=VALE9PT, E,POIS,oOD82, THB2,... 
Ro2eirTe2, NTBS2, NLUB2, NL982, NV82) 





vV8B3,A8B3,LE8B3,D8B3,LA8B3, BCS83=VALE(PI,POIS,0D83,TH83,... 
L83, NTT8 3, NTBO 3, NLL 3, NL98 3, NV8 3) 


BOI=1.0.B0 
BCII=2.0::(1.0+POIS)/E 
BI1=BOI+BC1I 
BI2=BI1 
B15=BI1l 
BI6=BI1 
BI7=BI1 
BI8=BI1 
BIJ=BI1 
BI10=BI1l 
BI11=BI1 
BI12=BI1 
BCDV8=BilsVDV8 
BCDV7=BI1*VDV7 
BCBV8=BI1#VBV8 
BCBV7=BI1+#VBV7 


BIS57=BOI+BCS5/V57 
BI68=BOI+BCS6/Vv68 
BI79=BOI+(BCS7+BCS71+BCS72+BCS7 3+BCS7H+BCDV7+BCBV7 )/V79 
BI8l= BOI+(BCS8+BCS81+B0S82+BCS6 3+BCSEU+BCDVS+BCBVS ) / Vs ~1 


CQ57,CQ57S=LINCO(A5,D5, D5) 
CQ58, CQ68S=LINCO(A6,D6,0D6) 
CQ79,CQ79S=LINCO(A7,D7,0D7) 
CO81, CQ81S=LINCO(A8, D8,0D8) 
@471,C011S=LINCO(A11,D11,0DR) 
CQ12,CQ12S=LINCO(A12,D12,0DR) 


P1Z=(PIN-PR)/2,0 
P2Z=012 
P57Z=P1Z 
P68Z=P1Z 
ey 9Z=P1Z 
POlZ=P1Z 
P117=P1Z 
P12Z=P1Z 


Q7Z=0.0 

@67=0.0 

@97=0 <0 

Q10Z=0.0 
Q11Z=0.0 
Q12Z=0.0 
Q57Z=0.0 
Q68Z=0.0 
Q79Z=0.0 
Q81Z=0.0 


Al1=PI*(D11+D11-D11R*«D11R)/4.0 
Al2=PI*(D12s+D12-D1l2R+D12R)/h.0 








D11=DLR 
D1l2=DLR 
D11R=DRR 
D12R=DRR 


ASS=A5::A5 
A6S=A6A6 
A7S=A7*A7 
AS8S=a8:-A8 


LS7=LA5 
A68=LA6 
L79=LA7 
L81=La8 


LEP57=LEP5 
LEP7 9=LEP7+LDV7+LBV7 
LEP68=LEP6 
LEP81=LEP8+LDV8+LBV8 


V57=VPS 
V68=vVP6 
V79=VP7+V7B1l+V7 B2+V7B3+VDV7+VBV7 
VO1=VP6+VOB1+V5B2+VEB3+VDVS+VBVE 


hpV7= SPECV (PDR7,QDR7; DDV7) 
LDV8=LDV7 
LBV7=SPECV( PBR7, RER7, DBV7) 
LBV8=LBV7 
VBV8=VBV7 
VDV8=VDV7 


DRE=ODRE-2.0+THRE 
BIR=BOI+2.0::((1.+POIS)*+ODREsODRE+(1. geile ae (DRE)/.<- 
( Es: (ODRE*ODRE-ORE*ORE) ) 


CEL=CELCC /( PCEL#98 3.1) 
CIL=CILCC/ (PCIL#98 3.1L) 


XRZ=0 .0 
ZRDZ=0 .0 


Dl=3.°-DLV 
D2=Dl 


LC=0.53:(LCY=-TR) 


DYNAMIC 
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#% PROCEDURE SERVO-VALVE ++ 


PROCED Q1,Q2,Q93;, QU=SERVO(AG, KeO,LAMtS. PIN, PR,P1,Pe,..- 
XV ,XVM,RHO1, RHO2, RHGR) 

CEX=KCO+EXP ( -LAMB:: (Xv/ XVM) #2) 
IF(XV)10, 20, 30 

10 Q1=CEX#(PIN-Pl) 
Q2=AG“SQRT (XVeXV%(PIN-P2)/RHO2) - 
Q3=AG#SQRT (ZV#ZV=(P1=PR) /RHOR) 
Qh=CEX*(P2=PR) - 
GO Tc oO 

20 Ql=CEX*(PIN-P1) 
Q2=CEX: ( PIN-P2) 
Q3=CEX* (P1l-PR) 
Qu=CEX+(P2-PR) 
GO TO Oo 

30 Q1=AGSXV#SQRT((PIN-P1)/RHO1) 
Q2=CEX:: ( PIN-P2) 
Q 3=CEX:: (P1-PR) 
Qu=AG(XV(SQRT( (P2-PR)/RHOR) 

ho CONTINUE 

ENDPRO 


“st PROCEDURE KPRO  +Hx% 


PROCEED K5,K6,K7,K8,K9,K1O=KPRO(P1,P2,P5,P6,P7,P8,P9,P10,K5P,... 
K5M,K6P,K6M,K7P,K7M,K8P,K8M,K9P,K9M,K1OP,K10M) 

IF(P1.GE.P5)GO TO 100 
K5=K5M 
GO TO 101 

100 K5S=K5P 

101 IF(P2.GE.P6)GO TO 102 
K6=K6M 
GO TO 103 

102: K6=K6P 

103 IF (P5.GE.P7)GO TO 104 
K7=K7M 
GO TO 105 

10h K7=K7P 

105 IF(P6.GE.P8)GO TO 106 
K=K8M 
CO. ro! 107 

106 K8=K8P 

107 IF(P7.GE. P9)GO TO 108 
K9=K9OM 
GO TO 109 

108 K9=K9P 

109 IF(P8.GE.P10)GO TO 110 
K10=K10M — 
GO TO 111. 

110 K10=K10P 

111 CONTINUE 

ENDPRO 





#% PROCEDURE EFFECTIVE LINE LENGTH = 


PROCED LES7=PRO57 (LEP57,L57,Q57, a RHOS7 ,MUS57) 
CLE57=@57:: -RHO57/ (MU57:: -DS) 
IF(CLE57.LT.785.l.)G0 TO 200 
LES7=L57+LEP57 ~ 


GO TO 201 
200 LES7=L57+1.27 3E-0 3% ‘CLES7# LEPST 
201 CONTINUE 
ENDPRO 


PROCED LE68=PR068 (LEP68, L68,Q68, DORHO68, MU68) 
CLE68=Q68RH068 / (MU68::D6 ) 
IF(CLE68.LT.785.4)GO TO 202 
LE68=L68+LEP68 


GO TO 203 
202 LE68=L68+1. 27 3E-0 3+ CLE68#LEP68 
203 CONTINUE 
ENDPRO 


PROCED LE79=PRO79(LEP79,L79,Q79,D7,RHO79,MU79) 
CLET7 9=Q7 9:+BHO7 9/ (MU79:+D7 ) 
IF(CLE79.LT.785.k4)G0O To 20h 
LE79=L79+LEP79 


GO TO 205 
20h LE7 9=L7 941. 27 3E~0 3%*CLE7 9::LEP79 
205 CONTINUE 
ENDPRO 


PROCED LE81=PRO081(LEP81,L81,Q81,D8, RHO81,MU81) 
CLE81=981+:RHO81/ (MU81+D8 ) 
MeCLEOl. LTo5.4)GO TO 206 
LES1=L81+LEP81 


GO TO 20/7 
206 LE81=L814+1. 27 3E-0 3%+CLE81*LEP81 
207 CONTINUE 
ENDPRO 


“x PROCEDURE MASS OF OIL «Hc 


PROCED MDOIL,MT=MOPRO(MR,XRD,RHO11,RHO12,V11,V12,V11D, 
vl2D, REO57 , RHO68,RHO79, REOS1, RHO1,RHO2,V1,V2 ,V57, ve, 
V79;5 V8l,... 
RHOR, VR) 
IF(XRD.GE.0.)GO TO 00 
MOIL=RHO11::V11+RHO7 92:V7 9+RHO5 7%: V57+RHO-V1+RHOR::VR 


GO To hOl 
4.00 MOIL=RHO12*V12+RHO81:-V81+RHO68+RHO2:: ‘V2+RHOR=VR 
MDOIL=RHO12*V12D 
hol MT=MR+MOIL 
ENDPRO 
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P1D=(Q1-Q3-Q5)/V1*B1I) 
P2D=(Q2-Qu.-Q5)/(VexBeT) 
P57D=(Q5-Q7) /(B571#V57) 
P68D=(Q6-Q8)/ ( B68I#V58) 
P79D=(Q7-Q9)/(B79I#V79) 
P81D=(Q8-Ql10)/(B81IVv81) 
P11p=(Q11-QIL-V11D)/(V11I«V11) 
Pl2D=(Q12+QIL+V1l2D)/(Bl2I#Vv1l2) 


PS=2.0%P57=-P7 
P7=2.0%P79-P9 
P9=(Q9-Q11)/CEL 
p6=2.0+p68-ps 
p8=2.0:P81-P10 
P10=(Q10-912)/CEL 
s+ LINE AND RAM FLOWS **< 


Q5=A5+SQRT (2.0% (P1=-P5)/(RHO5S#KS) ) 
Q6=A6:-SQRT(2.0%(P2-P6)/(RHO6*K6) ) 


THE FOLLOWING EQUATIONS IN IMPLICIT FORM WERE INSERTED IN 
AN EFFORT TO CLEAR THE ALGEBRAIC LOOP PROBLEM. IT DID NOT 
WORK. 


THE NEXT LOGICAL STEP WOULD BE TO COMPACT ALL EQUATIONS 
SUCH THAT EACH VARIABLE IS EQUATED ONLY TO ITSELF OR A 
VARIABLE WHICH IS DIRECTLY COUPLED TO AN INTEGRAL. 


Q7=IMPL(Q7Z,0.001,Q7CC) 

Q7C=(( P5-P7-CQ57S*RHO57#ABS (Q57 (-0.5%RHO7#K7Q707/Q5S) .-. 
/ (CQ57*#MUS7*LE57 ) )+2.-Q5 

Q7CCc=Q7C 

Q8=IMPL(Q8Z,0.001, Q8cc) 

Q8c=(( P6-P8-coQ68S:-RH068:-968+ABS (Q68)-0.5+RHO8:K8:Q8/A6S)... 
/ (6Q68++MU68:-LE68 ) )#2.-Q6 

Q8cc=98c 


Q9=IMPL(Q9Z,0.001,99CC) 

Q9C=( (P7-P9=CQ7 IS:RHO7 9+Q7 IHABS (Q79) -0. SH RHOI#*KIZQI“QI/ATS... 

» £(CQ79#MU79LE79) )#2.-Q7 

Q9CC=Q9C 

Q10=IMPL(Q102Z,0.001,Q10CC) 

Q10c=( ( P8-P10-cQ81S+RHO81::Q81+ABS (Q81)-0.5%RHO1O:+K10%+Q10+Q10... 
/A8S)/(CQ812MU81-:LE81 ) )+2.-Q8 

Q10cc=Q10¢c 


Q1l1=IMPL(Q11Z,0.001,Q11C) 
Ql1ic=(P9-P11-CQ11SsRHO11+Q11:-ABS (Q11) )/(CQ11#MU11:-L11 ) 
Q12=IMPL(Q127,0.001,012C) 
Q12c=(P10-P12-CQ125:+RHO1L2+Q12xABS (Q12) )/ (CQL2sMU12L12) 
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QS7=IMPL(Q572,0.0001,Q57CG) 

BS oe SeeRHOS 12051 #ABS (ST =O. SHRHOTHRTEQT#RT/A2S) + 
7(CQST#MUST# LEST ) 

-Q57C0C=Q57C 


Q68=IMPL(Q68Z,0.001,Q68CC) 

Soren CQbOSHRHOOOABS (Q6B)-0 - 5#RHOBHKEHQSEQE/ACS) «+ 
7( CQ68MU68%LESS ) 

Q68cc=Q68C 


Q79=IMPL(Q79Z,0.001,9790C) 

Bae oro a7 SSeRHOTIZQTIRABS (Q79)~0 + S#RHOIK IAI AI/ATS) +++ 
/ (CQ7 9%*MUT 9#LET 9) : 

Q79CC=Q79C 


Q81=IMPL(Q81Z,0.001,Q51CC) 

Mee (pb P10-CQBISHRHOB1#ABS (981) -0. 5#RHOLORKIOFALONONN. + 
7K88) / (CQ81#MU81::LES1 ) 

Q@81icc=Q81c 


QIL=CIL#(P11-Pl2) 


BLI=CFA/P1+BIL 
B2I=CFA/P2+Bie 
BRI=CFA/PR+BIR 
BOI=CFA/P5+B15 
BOI=CFA/PO6+BI16 
B7I=CFA/P7+BIL7 
BOI=CFA/P8+B18 
BIL=CFA/PI+BII 
BLOI=CFA/P1LO+BI10 
BiLI=CFA/PL1+BI1l 
B12I=CFA/P12+Blle 
BS7 1=CFA/P57+Bi57 
R6SI=CFA/PO8+BI16 8 
B791=CFA/P79+ BI7T9 
ROLI=CFA/PS1+B181 


sseeeDENSITY VALUES *t4* 


RHOR=8.«.256E-0548 .23E-05% PR*BRI 
ROH1=8.250E-05+8.423E-O5*P1sBlt 
RHOL=8. 256E- 0548.2 38-05% PeHECI 
RHOS=8 . 256E-05+8 1.2 3E-05#Po#BOL 
RHOG=B . 256E-05+8.423E-05%POXBOL 
RHO7=8.. 256E-05+8 .423E-O5"P7#BTI 
RHOB=8 .256E-05+8. 42 3E-05% POuBOL 
RHO9=8. 2568-0548. 23E-05% POXBII 
RHOLO=8 .256E-0548.423E-O5sPLOxBIOL 
RHO11=8 « 256E-0548 «2 3E-05#PLIHE III 
RHOL2=8 . 256E-05+8 2 3B-O5ePloxBlel 
RHOS7=8 . 250E-05+8 42 3E-OS#POTHBSTI 
RHOOB=B . 256-0548 .i2 3E-05# P68 BOOT 
RHO79=8 «2508-0548. 2 3E-O52P7 HBIIT 
RHOBI=8 , 250E-05+8.Y23E-O5*PS1ZBS1T 


OOO ee 





wit VISCOSITY EQUATIONS +x 


MUS7=8 .198E-06EXxP (P57+B571) 
MU68=8 .198E-06*EXP ( P68+B68T ) 
MU79=8.198E-06::EXP ( P79B791) 
Mu81=8 .198E-06EXP( P81+B81T) 
MU11=8 .198E-06+EXP(P11+B11I) 
MU12=8 .198E-06EXP(P12+B1l2I) 
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